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Using extensive Monte Carlo simulations we study the phase diagram of a symmetric binary 
(AB) polymer blend confined into a thin film as a function of the film thickness D. The monomer- 
wall interactions are short ranged and antisymmetric, i.e, the left wall attracts the A-component 
of the mixture with the same strength as the right wall the B-component, and give rise to a first 
. order wetting transition in a semi-infinite geometry. The phase diagram and the crossover between 

' different critical behaviors is explored. For large film thicknesses we find a first order interface 

localisation/delocalisation transition and the phase diagram comprises two critical points, which are 
the finite film width analogies of the prewetting critical point. Using finite size scaling techniques 
we locate these critical points and present evidence of 2D Ising critical behavior. When we reduce 
q , the film width the two critical points approach the symmetry axis cj> — 1/2 of the phase diagram 

Sand for D w 2R g we encounter a tricritical point. For even smaller film thickness the interface 
localisation/delocalisation transition is second order and we find a single critical point at <j> = 1/2. 

Measuring the probability distribution of the interface position we determine the effective inter- 
action between the wall and the interface. This effective interface potential depends on the lateral 
system size even away from the critical points. Its system size dependence stems from the large but 
finite correlation length of capillary waves. This finding gives direct evidence for a renormalization 
of the interface potential by capillary waves in the framework of a microscopic model. 



I. INTRODUCTION. 
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Confining a binary mixture one can profoundly alter its miscibility behavior, jl] ^) If a mixture is confined into a 
quasi one-dimensional (e.g., cylindrical) pore no true phase transition occurs, unlike the prediction of the mean field 
7-H , theory. Fluctuations destroy long-range order and only a pronounced maximum of the susceptibility remains in the 
^ ■ vicinity of the unmixing transition in the bulk. In a two-dimensional system (e.g., a slit-like pore or a film) with 
identical surfaces a true phase transition occurs (capillary condensation) and the shift of the critical point away from 
its bulk value has been much investigated. [0 The confinement changes the universality class of the transition from 
3D Ising critical behavior in the bulk to 2D Ising critical behavior in the film. The latter manifests itself in much 
flatter binodals in a film close to the unmixing transition than in the bulk. No such change of the critical exponents 
' is observed in mean field theory. 

The phase behavior of symmetric mixtures in a thin film with antisymmetric surface interactions has attracted 
"j***; i abiding interest recently. ]?|-^2]| The right surface attracts one species with exactly the same strength as the opposite 
surface attracts the other species. In contrast to capillary condensation, the phase transition does not occur close 
to the unmixing transition in the bulk, but rather in the vicinity of the wetting transition. Close to the unmixing 
transition in the bulk, enrichment layers at the surfaces are gradually built up and an interface is stabilized in the 
middle of the film. In this "soft-mode" phase the system is laterally homogenous - no spontaneous breaking of the 
symmetry occurs. If the wetting transition of the semi-infinite system is of second order one encounters a second order 
localisation-delocalisation transition slightly below the wetting transition temperature. The system phase separates 
laterally into regions where the interface is located close to one surface (localized state). The order parameter, i.e., 
. £^ | the distance between the interface and the center of the film, grows continously. This prediction of phenomenological 
theories has been corroborated by detailed simulation studies |l(],[l^,[l4|] and it is also in accord with experimental 
H I findings. |1|J| 

If the wetting transition is of first order and the thickness of the film not too small, mean field calculations p7 18 
predict the occurance of two critical points which correspond to the prewetting critical point of the semi-infinite 
system. Unlike the wetting transition, || the prewetting transition can produce a critical (singular) behavior in a 
thin film, because only the lateral correlation length diverges at the prewetting critical point; the thickness of the 
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enrichment layers at the surfaces remains finite. The mean field treatment invokes approximations and it cannot be 
expected to capture the subtle interplay between 2D Ising fluctuations at the critical points, "bulk" -like composition 
fluctuations, and interface fluctuations typical for the wetting transition. |ll| Consequently, a detailed test of the mean 
field predictions via Monte Carlo simulations is certainly warranted and elucidates the role of fluctuations. Using 
Monte Carlo simulations of the Ising model Ferrenberg et al. studied the interface localisation-delocalisation transition 
also for the case that the wetting transition of the semi-infinite system is of first order. Jl9| The simulation study 
was centered on the dependence on the film thickness, which is a convenient parameter to be varied in experiments. 
However, the study was restricted to the coexistence between strictly symmetric phases and many questions remained 
open. 

The general features of the phase behavior are shared by all binary mixtures. Here, we present large scale Monte 
Carlo simulations aiming at investigating the phase behavior of a symmetric binary polymer blend confined between 
antisymmetric walls. Computationally, simulations of a polymer blend [pof are much more demanding than studying 
simple fluids (e.g., the Ising model), but recent mean field calculations made detailed predictions for the phase behavior 
of confined polymer mixtures jl7],[l8) and serve as guidance for choosing the model parameters in the simulations. 
Simulating polymer blends, we can, at least in principle, control the importance of fluctuations by varying the degree 
of interdigitation, i.e., the chain length. [ p8|j2C| ] The mean field theory is expected to become accurate in the limit 
of infinite interdigitation. In a binary polymer blend the wetting transition occurs at much lower temperatures than 
the critical temperature of the unmixing transition in bulk. pl[ Hence, "bulk" -like composition fluctuations are not 
important in the vicinity of the wetting transition temperature and we can isolate the effect of interface fluctuation. 
Moreover, these systems are also suitable candidates to examine the phase behavior experimentally. Indeed, one of 
the first studies of the "soft-mode" phase has employed a binary polymer blend. [pUf 

Our paper is broadly arranged as follows: First, we present a phenomenological description of the phase behavior in 
a film with antisymmetric short range surface interactions. Using a standard model for the effective interface potential 
we calculate the phase behavior in mean field approximation, discuss the regime of validity of the mean field approach, 
and consider the crossover between the different critical behaviors. Second, we briefly describe our coarse grained 
lattice model for a binary polymer mixture. Then, we present our Monte Carlo results: We obtain the phase diagram 
for film thicknesses ranging from D = 1.1 R g to 7R g , where R g denotes the radius of gyration of the polymer chains, 
investigate the critical behavior and present evidence that interface fluctuations renormalize the effective interface 
potential. We close with a comparison of the phase diagram to the behavior of the bulk and of films with symmetric 
boundary conditions. 



II. BACKGROUND. 



Rather than describing the configuration of the system by the detailed composition profile across the film, much 
qualitative insight into the thermodynamics can be deduced from the effective interface potential. Below the bulk- 
critical temperature enrichment layers of the prefered components form at the surfaces and stabilize an AB interface 
which runs parallel to the walls. The effective interface potential g wa \\{l) describes the free energy per unit area as 
a function of the distance I between this AB interface and a wall. In the case of short range forces between the 
monomers and the walls, the interface profile becomes distorted in the vicinity of the walls and this gives rise to an 
interaction which decays exponentially as a function of the distance 1 between the AB interface and a single wall: 

<?wail(0 = aexp(— XI) — 6exp(— 2X1) + cexp(— 3X1) (1) 

This expression retains only the lowest powers of exp(— XI), which are necessary to bring about a first order wetting 
transition of the semi-infinite system. The coefficient a is explicitely temperature dependent, while the temperature 
dependence of b and c is neglected, c > is assumed throughout the discussion. All coefficients are of the same 
magnitude as the interfacial tension a between the coexisting bulk phases. For polymer blends this quantity scales with 
chain length N and monomer number density p like vN /Rg- N = (pR g /N) 2 measures the degree of interdigitation. 
1/A denotes the spatial range of the interactions and it is of the order R g . b < gives rise to a second order wetting 
transition at a — 0; and b — to a tricritical transition. For b > one encounters a first order wetting transition 
at a wct = b 2 /4c where the thickness of the enrichment layer jumps discontinuously from Z_ — l/Aln(2c/6) to a 
macroscopic value. Jll| The wetting spinodals take the values a > (from the wet phase) and a < b 2 /3c (from the 
non-wet phase). The concomitant prewetting line terminates at the prewetting critical point a pwc = 16a wc t/9 and 
Z pwc = l/Aln(9c/26). 

We approximate the effective interface potential in a film to be the linear superposition of the interactions originating 
at each wall and analyze the behavior. SCF calculations |l8| lend support to this approximation. The interface 
potential in a film of thickness D takes the form: 
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g(l) = 5waii(0 + 9wa,n{D - I) - 2.g wan (D/2) 

= 2aexp(-AD/2) (cosh (X[l - D/2]) - 1) - 26exp(-AD) (cosh (2X[l - D/2]) - 1) 

+2ccxp(-3AD/2) (cosh (3A[Z - D/2]) - 1) (2) 

In general, the phase boundaries depend on the variables a/c, b/c and XD. If we proceeded as in the Ref [fl3| by 
expanding the cosh in powers of [I — D/2], the further analysis would be rather cumbersome. A more transparent 
procedure employs the variable 

rh 2 = 2 exp(-AD/2) (cosh (X[l - D/2]) - 1) = (exp(-AD/4)A[Z - D/2]) 2 + higher orders of [I - D/2] (3) 

to rewrite the interface potential in the form 

m r~2/~2 n2 , .~2i - + u b-6cexp(-AD/2) a - a wct - frexp(-AD/2) 
gr(tj = c [m (m — r) + tm J with 7- = — and i = (4) 

The qualitative form of the effective interface potential has been inferred previously on the basis of a Landau expansion, 
Ji"8| here it is derived explicitely from the standard form of the interface potential ([[]) for a first order wetting transition 
in the semi-infinite system. Negative values of r correspond to second order localization-delocalization transitions, 
r = to a tricritical one, and positive values of r give rise to first order transitions, t measures the distance from the 
tricritical transition temperature (for r < 0) and denotes the triple temperature in the case of a first order interface 
localization-delocalization transition (cf . below) . For r < the phase boundaries depend only on the two parameter 
combinations r and t. In these variables the limit XD — > oo is particularly transparent: cr — ► 6/2, ct — ► a — a wot and 
rh — » exp(— XI). 



A. r < 0: second order and tricritical interface localization-delocalization transition 

A second order interface localization-delocalization transition (i.e. r < 0) will occur either if the wetting transition 
is second order (i.e., b < 0) or if the wetting transition is first order but the film thickness D small enough to comply 
with < b < 6cexp(— AD/2). This behavior is in accord with previous findings [||l9]|[|] and we shall corroborate 




this further by our present simulations. Since the coexisting phases are symmetric with respect to exchanging I and 
D — I, phase coexistence occurs at A/i coex = or dg/dl = (dg/drh)(dm/dl) — 0. From this condition we obtain for 
the binodals: 

Ai/4|r| for At < r 2 (2DMF) 

y^Kt/3 for At > r 2 (2DTMF) ^ 

The critical temperature is given by t c = —r 2 and At = t c — t denotes the distance from the critical temperature at 
fixed r. For r < the binodals at the critical point open with the mean field exponent /32DMF = 1/2. This corresponds 
to mean field critical behavior (2DMF) of a system with a single scalar order parameter, i.e., m = [//D — 1/2]. At larger 
distance the order parameter grows like m ~ (Ai)' 32D ™ F with /32DTMF = 1/4- The latter exponent is characteristic 
for the mean field behavior at a tricritical point (2DTMF). The crossover between mean field critical and tricritical 
behavior occurs around |Ai cross | ~ r 2 . As we decrease the magnitude of r — > we approach the tricritical point and 
the regime where mean field critical behavior is observable shrinks. At the tricritical point only the tricritical regime 
(2DTMF) exists, i.e., At cross = 0, and the binodals take the particularly simple form rh — (Arj/3) 1 / 4 . The crossover 
in the binodals for r = —0.4 is illustrated in the inset of Fig.|l|(a). 

Of course, the above considerations neglect fluctuations and the behavior close to the transition is governed by 
Ising critical exponents and two dimensional tricritical exponents, respectively. The crossover between Ising critical 
behavior (2DI) and tricritical behavior (2DT) occurs at |At cross | ~ j- 1 /^^^ where the crossover critical exponent is 
not 1/2 (as for the crossover between the mean field regimes) but rather 4/9. &|§ Following Ref we calculate 
the critical amplitudes and estimate the location of the crossover between mean field critical behavior and the region 
where fluctuations dominate the qualitative behavior. For small values of the order parameter m = [l/D — 1/2] we 
approximate rn w rncxp(AD/4)/(AD) and obtain for the mean field critical amplitudes: 

- exp(AD/4) - exp(AD/4) 

D 2DMF - - ^ and B 2DTMF = _ T7i -_ (6) 



The susceptibility of the order parameter above the critical temperature is related to the inverse curvature of the 
interface potential in the middle of the film 1/xD 2 = %^^ D / 2 ~ Using Eq(^) we obtain for critical and tricritical 
mean field transitions: 
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X = 2c(AD) 2 e M^D/2)M- 1 and G+ F = ~ 2 exp(AD/2) 7 mf = 1 (7) 

The ratio G^p/G^p °f the critical amplitudes above and below the critical point is universal and takes the mean 
field value 2 at the critical point and 4 at the tricritical point. At the transition the correlation length £m diverges. 
This lateral length is associated with fluctuations of the local interface position, i.e., capillary waves. In mean field 
approximation the parallel correlation length takes the form: 

-^fj = V^X hence |+ F = exp(AD/4) 7 mf = 1/2 (8) 

and £m F /£mf = anc ^ ^' res P ec tive-ly. 

Knowing the critical amplitudes we can estimate the importance of fluctuations via the Ginzburg criterium: |25| l 
As it is well known, mean field theory is self-consistent if the fluctuations of the order parameter in a volume of 
linear dimension £y are small in comparison to the mean value of the order parameter. For our quasi-two-dimensional 
system (d — 2) we obtain: 

x ! 2 fc 1 - 2 / d \ 2 \ d/2 , , ' f ^A^ 4 -^/ 2 for r<0 second order 



V °" / I A^ 3 ^)/ 2 for r=0 tricritical 

This result is as expected: For our quasi-to-dimensional system we obtain for a second order interface localization- 
delocalization transition At G« 2 di ~ | T"| exp(— XD/2)/Vn in accord with Ref [^3[ , while we obtain At <C Gz 2 dt ~ 
exp(— XD)/N upon approaching the tricritical point. For bulk (d = 3) tricritical phenomena Landau theory is 
marginally correct. 

Combining the above results we find the following behavior upon approaching the critical temperature: Far away 
from the tricritical point, i.e., r ^> exp(— \D/2)/V~N we find mean field tricritical behavior (2DTMF) for At 3> r 2 , 
mean field critical behavior (2DMD) for r 2 « Ai « \r\ exp(— AD/2), and finally two dimensional Ising critical 
behavior (2DI) for \r\ exp(— AD/2) ^> At. Closer to the tricritical point, i.e. r <C exp(— XD/2)/Vn, we find mean field 
tricritical behavior (2DTMF) for At > exp(-AD/4), two dimensional tricritical behavior (2DT) for exp(-AD/4) < 
At <C Cr 1 ^"""", and Ising critical behavior (2DI) for CY 1 /*"™ 3 ^ At. The prefactor G must be chose such that 
all crossover lines ( 2DI ^ 2DT, 2DT <-> 2DTMF, 2DTMF 2DMF, and 2DMF <-> 2DI) intersect in a common 
point. This yields C — (N exp(AD))~ 1+1 / 2< ^ croB3 . Of course, the term "crossover line" is not meant as a sharp division 
between different behaviors, but should be understood rather as a center of a smooth crossover region. Likewise, the 
above constant G may involve a constant of order unity which has been suppressed for simplicity. The two different 
sequences can be clearly distinguished in the Monte Carlo simulations, because the probability distribution of the 
order parameter exhibits a three peak structure |3(| only close to the tricritical point (2DT). We shall use this property 
to distinguish between the two different sequences in our MC simulations. The anticipated behavior is summarized 
in Fig.0(a). 

In the Monte Carlo simulation this rich crossover scenario is further complicated by finite size rounding. The Monte 
Carlo results are subjected to pronounced finite size effects whenever the correlation length becomes of the order of 
the lateral system size. In the mean field regime the correlation length scales like £|| ~ R g exp(AD/4)At -1 / 2 . Knowing 
the Ginzburg number for the crossover from 2DMF to 2DI behavior we estimate the correlation length in the Ising 
critical regime: Q 

£ o f\n/^A+-i/2iYA+/r<- \ f -R s exp(AD/4)Ai _1 / 2 for At > Gi 2D i , 1fV . 
£||~i? 3 exp(AD/4)A; /(Ai/G* 2DI ) - j ^| r |i/2#-i/4 At -i for Ai « Gi 2m (10) 

where we have assumed that the scaling function / assumes a power law behavior for small and large arguments and 
we have used the value f 2 Di = 1 appropriate for the divergence of the correlation length in the 2DI regime. Similarly, 
we determine the correlation length in the 2DT regime: 

t n m(\niA\\+-V*t(\+ir4 ^ _> I R a exp(AD/4)At- 1 /2 for At > Gi 2DT , in 

e„ ~ R g cxp(AD/4)At /(Ai/G* 2DT ) - j ^ exp(XD(3/4 _ ^^V^, At— for At « « 2DT (H) 



v%Ti = 5/9 denotes the exponent of the correlation length in the 2DT universality class. |22j-|2J] The correlation lengths 
at the various crossovers are compiled in Tab.l. The largest correlation length occurs at the crossover from 2DT to 
2DI behavior 
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^dt^di „ Rg exp(A ^ (3/4 _ ^ tri /20 cross ))7V 1 / 2 — i/ a *«~ | r |-^/*— (12) 

In order to observe the true Ising critical behavior for negative values of r, the system size L has to exceed this 
correlation length. In the vicinity of the tricritical point (i.e., for small negative values of r) this requirement is very 
difficult to be met in computer simulations. 



B. r > 0: first order interface localization— delocalization transition 

For positive values of r the interface potential exhibits a three valley structure. The three minima at m = i-s/r 
and to = have equal free energy at t = 0. This corresponds to the triple point. At lower temperatures an A-rich 
phase coexist with a i?-rich phase, and since the two phases are symmetrical the coexistence occurs at A/x coex = 0. 
The binodals below the triple point take the form: 




m = ±\— + \--- for t<0, r>0 (13) 



Above the triple temperature t > there are 2 two phase coexistence regions symmetrically located around to = 0. 
These phase coexistences terminate in two critical points. Since the coexisting phases correspond to a thick and 
a thin enrichment layer of the prefered phase at each wall, there is no symmetry between the coexisting phases, 
and the exchange potential A/i coox at coexistence differs from zero. Unfortunately, the phase boundaries for t > 
and r > depend not only on r and t but also on XD explicitely, and we have determined them numerically. The 
dependence of the critical temperature t c on r for several values of XD is presented in Fig.[l](b). The coexistence curve 
for b/c= AAA and various values of XD are presented in the inset of Fig||(b). As the film thickness is decreased the 
critical temperature decreases and the critical points move closer to the symmetry axis of the phase diagram. They 
are determined by the condition: 

d 2 g d 3 g 

~~dP = ~dP = ° at t = tc to = to c (14) 

In two limiting cases a simple behavior emerges: 

(i) If \X{1 — ZD/2) j C 1 we can replace the derivative w.r.t. I by derivatives w.r.t. m and obtain: ]l8| t c = 7r 2 /5 and 
to c = ±y/2r/5. This approximation holds for r -C exp(— AD/2). Expanding g in powers of dm — m — rh c we obtain 
(omitting an irrelevant term linear in dm): 

7c 

g(m) w c(t c — t) Sm 2 + Arc 8m A + — Sm 5 + c Sm e (15) 

o 

This allows us to calculate the binodals in the vicinity of the critical points, the susceptibility and parallel correlation 
length. The presence of a 5th order term in Sm in the expansion ( |T5| ) is a manifestation of the fact that the phase 
boundaries of the prewetting transitions are not symmetric around rh c . This lack of symmetry is also evident from the 
numerical results in Fig.Q(b) (insert). The critical amplitudes scale in the same way with r, N, and XD as for r < 0. 
In particular we find for the crossover between 2DMF behavior and 2DI behavior Gi^di ~ \r \ exp(— XD/2)/Vn. 

(ii) In the limit of large film thickness XD — > oo, the critical point tends towards the prewetting critical point at 
t c = tpwc = 7r 2 /9. In this limit confinement effects are negligible and the coexistence curves in the vicinity of the 
critical points corresponds to the prewetting lines at the corresponding surfaces. We expect the same critical behavior 
as at the prewetting critical point. In this case the Ginzburg number does not depend on the film thickness. For 
XD — > oo we employ the interface potential at a single wall, and find for the validity of the mean field description: 

— ~ Qpwc » — ~ N- 1 ' 2 ~ Gi 2D1 for r>0 XD - oo (16) 



III. THE BOND FLUCTUATION MODEL AND SIMULATION TECHNIQUE. 

Modeling polymeric composites from the chemical details of the macromolecular repeat units to the morphology 
of the phase separated blend within a single model is not feasible today even with state-of-the-art supercomputers. 
Yet, there is ample evidence that by a careful choice of simulation and analysis techniques, coarse grained models 
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of flexible polymers - like the bond fluctuation model |2(J,|2(J - provide useful insights into the universal polymeric 
features. In the framework of the bond fluctuation model each effective monomer blocks a cube of 8 neighboring 
sites from further occupancy on a simple cubic lattice in three dimensions. Effective monomers are connected by 
bond vectors of length 2, y/S, VE, 3, or y/lO in units of the lattice spacing. The bond vectors are chosen such that 
the excluded volume condition guarantees that chains do not cross during their motion. 27 Each effective bond 
represents a group of n « 3 — 5 subsequent C — C-bonds along the backbone of the chain. |28| Hence, the chain 
length N = 32 employed in the present simulations corresponds to a degree of polymerization of 100 — 150 in a real 
polymer. If we increased the chain length N, the mean field theories would yield a better description of the equilibrium 
thermodynamics (self-consistent field theory is believed to be quantitatively accurate in the limit N — > oo) but the 
length scale of the ordering phenomena would be larger. Hence, our choice of N is a compromise determined by the 
computational resources. The statistical segment length b in the relation for the radius of gyration R g — byj N/6 is 
b = 3.05 (i.e., R g « 7 for N = 32). 

We study thin films of geometry L x L x D, Periodic boundary conditions are applied in the two lateral directions, 
while there are hard impenetrable walls at z = and z = D + 1 modeling a film of thickness D. The average number 
density in the film is po = I/I6, i.e., half of the lattice sites are occupied by corners of monomers. This density 
corresponds to a melt or concentrated solution. The density profile of occupied lattice sites, normalized by the bulk 
value, is presented in Fig.|] for film thickness D = 24 and 48. For this choice of temperature and monomer-wall 
interaction an interface is stabilized in the center of the film. Due to the extended shape of the monomers and the 
compressibility of the fluid there are packing effects at the walls. |2lf| Overall the walls are repulsive and the monomer 
density is slightly reduced in the boundary region. The spatial extension of this region is independent of the film 
thickness. Moreover, the density is reduced at the center of the interface as to reduce the energetically unfavorable 



contacts between unlike species. [ 
cause the density in the "bulk"-li 



Both effects are not incorporated into the mean field calculations JL7 18 and 



te region of the film to be slightly larger for thinner films than for thicker ones. In 
the following we employ the density of occupied lattice sites in the layers 5 < z < 8 as a measure of the density of 
the film. For large D the data are compatible with a behavior of the form p = po(l + 0.85/ D). The film thickness 
ranges from D — 12 w 1.7 R g to D = 48 w 7R g and we vary the lateral extension over a wide range 48 < L < 264 to 
analyze finite size effects. In the two layers nearest to the walls, monomers experience a monomer-wall interaction. 
An A-monomer is attracted by the left wall and repelled by the right wall, the interaction between i?-monomers and 
the walls is exactly opposite. Each monomer-wall interaction changes the energy by an amount e w — 0.16 in units 
of UbT ' . For these parameters the wetting transition and the phase diagram of a blend confined between symmetric 
walls has been investigated previously. pl[ 

Binary interactions between monomers are catered for by a short ranged square well potential e = —£aa = —£bb = 
€ab = 1/kgT which is extended up to a distance %/6- The phase separation is brought about by the repulsion between 
the unlike species. The Flory-Huggins parameter is x — 2z fje where z e s ~ 2.65 denotes the effective coordination 
number in the bulk |30||2C| ] at pp = 1/16. For e w — 0.16 previous simulations find a strong first order wetting transition 
at T wct = 1/fiwct = 14.1(7). Q This value corresponds to \N ~ 12, which is well inside the strong segregation limit. 

The polymer conformations are updated via a combination of random monomer displacements and slithering snake- 
like movements. The latter relax the chain conformations about a factor of N faster than the local displacements. ]3C| ] 
We work in the semi-grandcanonical ensemble, |3l|| i.e. , we control the temperature T = 1/e and the exchange potential 
Ap, between the two species, and the concentration fluctuates. This semi-grandcanonical ensemble is realized in the 
Monte Carlo simulations via switching the polymer identity A ^ B at fixed chain conformation. The different Monte 
Carlo moves are applied in the ratio: slithering snake : local displacements : semi-grandcanonical identity switches = 
12:4:1. During production runs, we record all 150 slithering snake steps the composition, the energy, and the surface 
energy and obtained the joint probability distribution in form of a histogram. We use the semi-grandcanonical identity 
switches in junction with a reweighting scheme, ]29| , |32]] i.e., we add to the Hamiltonian of the system a reweighting 
function H lw = H OT i g +W \4>) , which depends only on the overall composition = n^ oly / (n^ oly +n® oly ). n^ oly and n^ oly 
denote the number of A and B polymers in the simulation cell, respectively. The choice W(<j>) « — lnP(0) , where 
P(4>) denotes the probability distribution of the composition in the semi-grandcanonical ensemble, encourages the 
system to explore configurations in which both phases coexist in the simulation cell. Otherwise these configurations 
would be severely suppressed due to the free energy cost of interfaces. In the framework of this reweighting scheme the 
system "tunnels" often from one phase to the other and this allows us to locate the phase coexistence accurately and 
measure the free energy of the mixture as a function of the composition cf>. Use of histogram extrapolation technique 
J33| permits histograms obtained at one set of model parameters to be reweighted to yield estimates appropriate to 
another set of model parameters. We employ this analysis technique to obtain estimates for the reweighting function 

w{4>). 
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IV. RESULTS. 



Firstly, we locate the critical points of the phase diagrams. For very small film thickness we find a second order 
localisation-delocalisation transition even though the wetting transition is of first order. Swift et al have predicted this 
behavior in the framework of a s quare g radient theory Q and such a behavior is also born out in our self-consistent field 
calculations for polymer blends ^Jlg ] and simulations of the Ising model. |nj Upon increasing the film thickness we 
encounter a nearly tricritical transition. A truly tricritical transition cannot be achieved by tuning the film thickness 
only, because of the discreteness of the lattice, but it could be brought about by varying the monomer-wall interaction. 
In an experiment using real materials, of course, the film thickness can be varied continously, and a truly tricritical 
transition is in principle always accessible. For even larger film thickness the interface localisation-delocalisation is 
first order and we find two critical points at <j> ^ 1/2. 

Secondly, we locate the triple line for the two largest values of the film thicknesses and discuss how capillary waves 
lead to a strong dependence of the effective interface potential on the lateral system size. 

Thirdly, we detail our results on the thickness dependence of the phase diagram and relate our findings to the 
binodals of the bulk and the mixture confined into a film with symmetric boundaries. 



A. Critical points. 



1. D = 8a l.lRg and D — 12 

For film thicknesses which are comparable to the radius of gyration of the molecules, the effective interface poten- 
tial originating from the two surfaces strongly interfere. This might change the order of the interface localisation - 
delocalisation transition from first to second. In this case, a single critical point occurs on the symmetry axis <f> = 1/2 
of the phase diagram. The transition is thought to belong to the 2D Ising universality class. In Fig.|| (a) we present 
the probability distribution of the composition for various inverse temperatures e, film thickness D = 8 and lateral film 
extension L = 80. Upon increasing the monomer-monomer interaction e the probability distribution P(</>) changes 
from single-peaked to bimodal, which indicates that a phase transition occurs in this temperature range. No signature 
of the trimodal distribution occurs and, hence, we conclude that the system is far away from the tricritical point, i.e., 
|r| > exp(— \D/2)/V~N. In this case, we expect a crossover from 2DMF to 2DI behavior. 

Along the coexistence curve Afi = and its extension to higher temperatures we use the cumulant intersection 
method to locate the critical point. Bq] In the vicinity of the critical point the probability distribution of the order 



parameter m = <f>— cj) coex = 0—1/2 scales to leading order like: 35 



P(m, L, t) ~ L p / v P*{L f} /' J m, L^t) (17) 

where t = (e c — e)/ e c denotes the distance from the critical point along the coexistence curve and (3 and v are the 
critical exponents of the order parameter and the correlation length. P* is characteristic of the universality class 
and has been obtained from simulations of the Ising model pif at the critical temperature t — 0. Cumulants of 
the form {m 2 )Um\) 2 are expected to exhibit a common intersection point for different system sizes L at the critical 
temperature. B5] The value of the cumulant at the intersection point is universal. Our simulation data are presented 
in panel (b) and exhibit some corrections to scaling due to the crossover 2DMF to 2DI behavior. Similar corrections 
were observed in simulations of a second order interface localization/delocalization transition in the Ising model. |l3| ] 
From the intersection points of neighboring system sizes and from the intersection of the cumulant with the universal 
value of the Ising model we estimate the critical temperature to be e c = 0.0520(5). 

In the inset of Fig.[|(b) we show the probability distribution normalized to unit variance and norm at our estimate of 
the critical temperature e c = 0.052, and compare the distribution to universal scaling curve of the 2D Ising universality 
class. The probability distributions for the smaller system sizes are slightly broader than the universal scaling curve, 
but the deviations decrease as we increase the system size. 

The simulation data for D = 12 are presented in Fig.||(c) and (d). As we lower the temperature the probability 
distribution of the composition for L = 48 changes from single peaked to bimodal. At intermediate values of e, however, 
a three-peak structure is clearly discernable. This is characteristic of the 2DT regime and indicates the vicinity of 
the tricritical point. In the phenomenological considerations this regime occurs only for \r\ < exp(— \D/2)/~\/n. We 
note that the distribution for that small lateral system sizes resembles at no value of e the universal shape of the 
order parameter distribution of the 2D Ising model. We conclude that the finite size rounding for this lateral system 
size sets in before we observe the crossover from 2DT to 2DI behavior, i.e., the correlation length ^ DT ^ 2DI m Eq 

L2T) exceeds the lateral systems size L. For such small lateral extensions the universal properties of the transition 
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are completely masked. Larger system sizes and a careful finite size scaling analysis is indispensable to determine the 
type of transition and accurately locate the transition temperature. 

The temperature dependence of the cumulant is presented in Fig.^|(d). There is no unique intersection point and 
the value of the cumulants at the crossing is larger than the universal value of the cumulant of the Ising class. This 
behavior indicates pronounced corrections to scaling due to the crossover from 2DT behavior away from the critical 
point to 2DI behavior at the critical point. From the intersection points of neighboring system sizes and from the 
intersection of the cumulant with the universal value of the 2D Ising model we estimate the critical temperature to 
be e c = 0.0589(10). 

The inset of panel (d) compares the distribution of the order parameter at our estimate of the critical temperature 
and the Ising scaling function. As we increase the lateral system size the "third" peak in the distribution vanishes 
and P{4>) gradually approaches the universal scaling curve. This indicates that our largest system sizes exceed the 
correlation length at the crossover from 2DT to 2DI behavior. The comparison of P(<fi) with the universal scaling 
curve for several system sizes accurately locates the critical point and gives evidence that the transition belongs to 
the 2D Ising universality class. 

For D < 12 we find a single interface localisation-delocalisation transition of second order at (j> = 1/2. 



2. D = 14 » 27? 9 : tricritical interface localisation-delocalisation transition 

The three-peak structure in the probability distribution for D = 12 and small lateral extensions L has indicated 
the vicinity of the tricritical interface localisation-delocalisation transition. Increasing the film thickness we need 
larger and larger lateral extensions to observe the 2DI behavior as ^ DT ^ 2DI diverges. Right at the tricritical point 
the distribution of the composition is expected to exhibit a three-peak structure for all lateral system sizes and the 
distribution, when scaled to unit variance and norm, coincides with a universal scaling function. Wilding and Nielaba 
|fj6f have obtained this scaling function via simulations at the tricritical point of the spin-1 Blume-Capel model |37]] 
in two dimensions. Assuming that the tricritical interface localisation-delocalisation transition belongs to the same 
universality class, we vary the film thickness D and the interaction strength e as to match the probability distribution 
of the composition onto the predetermined scaling function of the tricritical universality class. This strategy largely 
facilitates the search of the tricritical interface localisation-delocalisation transition. Fig.||(a) displays the probability 
distribution of the composition for film thicknesses ranging from D — 12 to D — 18 and the universal scaling curve. 
The temperature was adjusted for each film thickness such that the relative heights of the central and outer peaks 
correspond to the ratio of the universal scaling curve. For small D < D tY i the "valley" between the peaks is too 
shallow and for D > Dm the probability between the peaks is too small. For D 3> Dt rl this situation corresponds to 
the triple point (cf. below) and the probability of finding a system between the peaks is suppressed by the free energy 
cost of the interface between the phases with composition close to or 1 and the "soft-mode" phase with composition 
4> = 1/2. As we increase the film thickness the temperature at which the ratio between the peak height equals 1.2 
shifts towards lower temperatures and approaches the wetting transition temperature from above. 

Panel (a) of Fig.|] suggests that the tricritical transition occurs close to the film thickness D — 14. This is further 
corroborated in FigH(b), where we show the distribution function at e = 0.06151 for various system sizes. Within the 
statistical accuracy of our data the distribution functions for the larger systems sizes collapse well onto the universal 
scaling curve. For smaller systems the outer peaks are slightly sharper and centered at smaller values of the order 
parameter. Of course, no perfect data collapse can be expected because we can tune the film thickness only in units 
of the lattice spacing. In view of the statistical accuracy and possible systematic corrections to scaling, however, we 
did not attempt to vary the monomer-wall interaction e w as to achieve a better collapse. For D = 14 the system is 
very close to the tricritical transition. 



3. D = 24 « 3.5R g and D = 48 « 7R g : critical points for 4>^l/2 



Though the system is strictly symmetric the critical points for larger film thickness (D > Z?tri) do not occur at 
(j> = 1/2 but rather there are two critical points at critical compositions <j) c and 1 — (f> c . These critical points are the 
finite film thickness analogs of the prewetting critical points, which occur in the limit D — ► oo. ]l7 ] Below the critical 
temperature the phase diagram comprises two miscibility gaps. The coexisting phases correspond to surfaces with a 
thin and a thick enrichment layer of the preferred component. Due to the missing symmetry between the coexisting 
phases the coexistence value of the chemical potential A/i cocx differs from zero. We determine A/i coex via the equal 
weight rule, f$§] i.e., we adjust A/z such that 
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I f 1 f 1 

d<t> P{<t>) = d<f> P{4>) and <f>* = d<j> P(</>)0 (18) 
J<P* Jo 

Along this coexistence curve and its finite-size extension to higher temperatures we use the cumulant intersection 
to locate the critical temperature. This is shown in Fig.|5|(a) for the film thickness D = 24. For the system sizes 
accessible in the simulations the intersection points between cumulants of neighboring systems sizes systematically 
shift to lower temperatures and the value of the cumulant at the intersection point gradually approaches the value of 
the 2D Ising universality class from above. The latter is indicated in the figure by the horizontal line. From these 
data we estimate the critical parameters to be e c = 0.061(1), and 4> c — 0.18(2) and <fi c = 0.82(2) respectively. This 
corresponds to a critical thickness l c = D<f) c = 0.62R g of the enrichment layer. A similar procedure has been employed 
to locate the critical temperature in the film of thickness D = 48. The temperature and system size dependence of the 
cumulants are displayed in Fig.||(c). From this we extract the estimate e c = 0.0625(10) for the critical temperature 
and (f> c — 0.09(2) and <fi c = 0.91(2) for the critical compositions. This value corresponds to a distance between then 
wall and the interface of l c — 0.63R g . Since increasing the film thickness from 3.5R g to 7R g does not change T c or l c 
substantially, we are in the regime XD 3> 1 and the critical behavior is characteristic of the prewetting critical point 
in the semi-infinite system. 

The behavior of the cumulants and the very gradual approach of the probability distribution towards the Ising 
curve indicate pronounced corrections to scaling. For the simulation of the bulk phase diagram [fj9| a nice cumulant 
intersection has been obtained with system sizes in the range 24 3 to 56 3 . In the present study we employ systems 
with about an order of magnitude more polymers and obtain no clear intersection of the cumulants! There are three 
reasons for strong corrections to the leading 2D Ising scaling behavior: (i) The aspect ratio D/L of our simulation 
cell is always finite. Truly two-dimensional behavior can only be observed for vanishing aspect ratio, and our data 
might fall into the broad crossover region between three-dimensional critical behavior and the two-dimensional critical 
behavior. Such a crossover has been studied in our polymer model for neutral walls f40| and walls, which attract both 
the same species (i.e., capillary condensation). J2lj However, we note that unlike these situations there is no three- 
dimensional critical point in the vicinity for antisymmetric boundary conditions. The temperature of the unmixing 
transition in the bulk is a factor 4 higher than the critical point in a thin film. Since the critical point in a thin film is 
related to the prewetting transition of the semi-infinite system, i.e., a transition with no three dimensional analogy, we 
expect the corrections to be qualitatively different from the case of neutral or symmetric boundaries, (ii) Unlike the 
situation for small film thickness D = 12 the probability distribution of the order parameter is asymmetric, because 
the critical point does not lay on the symmetry axis of the phase diagram. This missing symmetry between the two 
phases gives rise to field-mixing effects, [|4j which manifest themselves in corrections of relative order £~( 1 ~ a ~/ 3 )/ zy '. 
These corrections are antisymmetric to leading order and, hence, are not expected to influence even moments (like 
the cumulants) of the order parameter distribution profoundly. The effects are, however, detectable in the order 
parameter distribution which we present in FigJ^(b) and (d). The distribution functions at our estimate of the 
critical temperature clearly lack symmetry and approach very gradually the symmetric scaling curve of the 2D Ising 
universality class, (iii) Additionally, there are corrections to scaling by non-singular background terms. One source 
of (non-critical) composition fluctuations are "bulk-like" fluctuations in the ^4-rich and £?-rich domains. In a bulk 
system, i.e., with periodic boundary conditions in all directions, the susceptibility is rather small. At e = 0.065 it takes 
the value x^ ulk = V (A<fi 2 ) = 0.047 with Acj) = </>—((/>). In a system of size 96 x 96 x 24 this susceptibility corresponds 
to composition fluctuations of the order \J (A<p 2 ) ~ 5 • 10~ 4 . Therefore, we believe that "bulk-like" composition 
fluctuations are not the major source of background terms. However, we cannot rule out that the presence of an AB 
interface gives rise to enhanced composition fluctuations. Another source of corrections to scaling stems from the 
fluctuations in the average interface position itself. Since the effective interaction between the interface and the wall 
is rather weak, they give rise to a finite but large susceptibility away from the critical point. We have estimated the 
susceptibility from the curvature of lnP(</>) close to the triple point (i.e., T m 0.9T C ), and we have obtained values 
of the order \t ~ 3 • 10 2 (and a smaller value is obtained if the interface is close to a wall.) For the same system 
size as above, this yields composition fluctuations of the order yj (Ac/) 2 ) ~ 0.04 (a value which should be compared 
to 4> C (D = 24) = 0.18(2)). This observation partially rationalizes why the peak in the probability distribution of the 
composition close to = 1/2 is always broader than the peak which corresponds to the phase in which the interface 
is close to the wall. As we approach the critical temperature composition fluctuations grow. At the critical point 
the typical composition fluctuations are of the order y Ac/) 2 ~ V Lpl v ~ d ~ i -1 / 8 , where we have used the critical 
exponents for the susceptibility 7 = 7/4 and the correlation length v = 1 appropriate for the 2D Ising universality 
class. Hence, for small system sizes typical fluctuations yield compositions which differ substantially from the critical 
composition; only for very large sizes the composition fluctuates in the vicinity of the critical value. Moreover, the 
critical density is much displaced from the symmetry axis 0=1/2 and typical fluctuations in a finite system are 
cut-off by the constraint < <j> or <j> < 1. Therefore, the susceptibility of a small system is reduced compared to the 
value expected from the leading scaling behavior. This observation is in accord with our Monte Carlo data, and a 
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similar reasoning has been used by Bruce and Wilding |41| in discussing background terms to the specific heat and 
the concomitant corrections to scaling in the energy distribution. 

B. The triple point. 

For the largest two film thicknesses D = 24 and D = 48 the interface localisation — delocalisation transition is first 
order and the concomitant two miscibility gaps join in a triple point. At this temperature an A-rich phase, a £?-rich 
phase and a phase where the interface is located in the middle of the film (4> = 1/2) coexist. The coexisting phases 
correspond to three peaks in the distribution of the composition. Upon increasing the lateral system size the peak 
positions do not shift (as opposed to the behavior at the tricritical point), the peaks become more pronounced and 
configurations with intermediate compositions are more and more suppressed, because of the presence of interfaces 
between the coexisting phases. 

The composition of the system and the average interface position are related via I = <pD (integral criterium) , where 
we assume that the coexisting bulk phases are almost pure, i.e., sa or 1. From the probability distribution we 

then calculate the effective interface potential g(l): 

9(l) = -^lnP(<i> = l/D) (19) 

In principle, not only fluctuations of the interface position (Al 2 ) but also "bulk" -like fluctuations contribute to 

bulk 

composition fluctuations (A<f> ) ~ jp{Al ) + jji D ■ Since the wetting transition in a binary polymer blend occurs far 
below the critical point of the bulk, the bulk susceptibility is very small, and the latter contribution can be neglected. 

The dependence of the free energy per unit area on the position of the interface is a key ingredient into the theory of 
wetting. 1 4^]43|,||,[l4] [i7) The interface interacts with the boundaries and the (bare) interface potential exhibits three 
minima. These correspond to the three coexisting phases. In the two phases with <p close to and 1, the interface 
is localized close to the wall, the interaction between the wall and the interface is rather strong, and the effective 
interface potential possesses a deep minimum. In the "soft-mode" phase the interface is only weakly bound to the 
center of the film and the minimum is much broader. In FigJ| we present the effective interface potentials for film 
thicknesses D — 24 (a) and D = 48 (b) and various lateral system sizes in the vicinity of the triple temperature. The 
three minima are clearly visible, however, the shape of the interface potential and the value of the minima depend on 
the lateral system size L. Moreover, the minima which correspond to the localized states broaden and (slightly) shift 
to larger distances between wall and interface upon increasing L (cf. inset). 

Fluctuations of the local interface position, i.e., capillary waves, lead to a renormalization of the effective interface 
potential g(l) and cause the dependence of g(l) on the lateral system size, which we observed in a microscopic model of 
a polymer mixture. Describing the configuration of the system only via the local position l(x, y) of the interface (sharp 
kink approximation) we write the coarse grained free energy in form of the capillary wave Hamiltonian: p| , ^6|j4^ ] 

H[l]=jd 2 x {|(VZ) 2 + 5 (0} (20) 

where a approaches the AB interface tension between the coexisting bulk phases for large separations between the 
wall and the interface. An increase of a at smaller distances I as revealed by previous MC simulations is neglected. 
pl[ In the vicinity of a minimum of g(l) we may approximate the interface potential by a parabola. 

g(Sl) = const + ^akjSl 2 (21) 

SI denotes the deviation of the local interface position from the position where the g(l) attains its minimum. £|| = 2w/k\\ 
is the parallel correlation length of interface fluctuations. For lateral distances much smaller than £|| the fluctuations 
of the local interface position are hardly perturbed by the interaction between the interface and the wall; the interface 
behaves like a free interface. For lateral distances which exceed £n capillary waves are strongly suppressed. £|| is larger 
for the minimum of g(l) in the center of the film than for the minima, in which the interface is localized at a wall. From 

the curvature of the effective interface potential g(l) for film thickness D = 24 we estimate k\ = J /aD 2 = 0.26 

and k% — 0.031, where we have used the bulk value a = 0.0382 for the interfacial tension at e = 0.068. For the thicker 
film we obtain k\ = 0.3, but the curvature in the middle of the film could not be accurately estimated. The value is 
of the order k^ ~ 0(0.005), and we expect this value to decrease exponentially with the film thickness. Hence, this 
fluctuation effect is the stronger the larger the film thickness. For the system sizes employed in the MC simulations 
kuL is of order unity. 
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In our Monte Carlo simulations the finite lateral system size L acts as an additional cut-off for the spectrum 
of interface fluctuations jTiJ and upon increasing L we extend the spectrum of interface fluctuations. Allowing for 
interface fluctuations we decrease the free energy of the system. Therefore, we expect the free energy density of the 
system to decrease when we increase the lateral system size, and we expect the effect to be the stronger the larger 
Therefore, the free energy of the "soft-mode" phase becomes smaller compared to the free energy of the phase, where 
the interface is located close to a wall when we increase L. This effect is clearly observed in the MC simulations. To 
be more quantitative, we consider a system where the laterally averaged interface position is at the minimum of g(l), 
and we expand the deviation Sl(x,y) from the minimum in a Fourier series 

oo 

dl(x, y) = ^ { anm cos ( ( ?n a; ) cos(q m y) + b nm cos(q n x) sm(q m y) + c nm sm(q n x) cos(q m y) + d nm sm(q n x) sin(q m y)} 

n,m— 

(22) 

with q n = 2nn/L. The coefficients aoo = &oo = c oo = <^oo — &o m = c 0m = d>o m = d n o vanish identically, all other 
coefficients can take any real value. Using this expansion (|2|) and the effective interface Hamiltonian (pp|), we calculate 
the average size of fluctuations 
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where the factor rj nm takes the values 7700 = 0, r/nQ = rjom = 1/2 and Tj nm = 1 for n / and m =/= in order to 
account for the restriction on the coefficients a, b, c, d. The additive constant is independent of the wavevector cut-off 
fcii . The dependence of the free energy on the system size is dominated by the small q behavior. In this regime the 
discrete nature of the wavevector space matters and, hence, we do not replace the sum over q by integrals. Using 
the measured values of the wavevector cut-offs we calculate the lateral system size dependence of the free energy 
difference between the "soft-mode" phase and the delocalized state. The results are compared to the MC data in 
Fig.^. Good agreement is found for large L, whereas there are deviations for smaller L. For small L the amplitude of 
the fluctuations becomes large and the a parabolic interface potential is no longer a good approximation - especially 
for the localized state where the interface is located very close to the walls. We have used histogram extrapolation to 
adjust the temperature such that the difference Ag(l) = g^ —Q\ of the minima vanishes. This corresponds to the equal 
height criterium for the triple point. The equal weight condition, which we have applied to determine the binodals 
close to the critical points, would require: Ag = -^j hi Both conditions agree, of course, when we extrapolate our 
results to L — ► 00. From this procedure we obtain the following estimates for the triple point: 1/etripic = 14.7(4) and 
triplc = 0.015, 0.5, 0.985 for D = 24 and l/e t = 14.2(4) and tr i P ic = 0.0066, 0.5, 0.9934 for D = 48. The thickness 
of the microscopic enrichment layer at the wetting transition temperature is of the order l wct = 0.05i? g ; a value which 
is consistent with expectation for strong first order wetting transitions. 

The dependence of the critical temperature and the triple temperature on the film thickness is summarized in Figj^. 
When we increase the film thickness the critical temperature l/e c shows a non-monotonic dependence. At D = 14 the 
tricritical point (where the critical temperature and the triple temperature merge) occurs at e t ri = 0.0615(5), for film 
thickness D = 24 we find e c = 0.0610(10) and at D = 48 e c = 0.0625(10). This effect is rooted in two opposing effects. 
On the one hand, the self-consistent field calculations predict the Flory-Huggins parameter Xc (D) to decrease 
upon increasing the film thickness D for an incompressible fluid. This shift in temperature decreases exponentially 
with the film thickness. One the other hand, packing effects, which are not incorporated in the self-consistent field 
calculations, increase the density in the "bulk" -like portion of the film when we decrease the film thickness. These 
packing effects at the walls depend strongly on the computational model, but qualitatively similar effects might occur 
in experimental systems as well. This thickness dependence of the density in the middle of the film modifies the 
relation between the depth of the square well potential and the ^-parameter. This leads to a behavior of the form 
e c ~ Xc CF /(l + 0. 85/13), where we use the dependence of the density profile (cf. Fig.||) on the film thickness as 
obtained by direct measurement in the Monte Carlo simulations. A dependence of the fluid packing structure on 
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the density is neglected. A similar 1/D correction to the difference in surface free energies between the A and the 
B-rich phase has been observed in previous simulations. |2lf| . Attempting to separate these two effects we also present 
[(1 + 0.85/ D)e c ]~ 1 which corresponds to the inverse Flory-Huggins parameter. Within the error bars the behavior of 
this quantity is consistent with the mean field prediction. The critical value of the inverse Flory-Huggins parameter 
increases and the triple value decreases as we increase the film thickness. The latter approaches the wetting transition 
temperature [^l) T wct = 14.1(7) from above. 

C. The phase diagram. 

For film thickness D = 48 we have determined the complete phase diagram. Close to the critical point we assume 
2DI behavior with an exponent (3 = 1/8 for the order parameter and employ finite size scaling to estimate the critical 
amplitude. Outside the critical region but above the triple temperature we have estimated the location of the binodals 
via the equal weight criterium in a system of size L — 64; but no finite size analysis has been applied. The phase 
diagram for a blend confined into a film with antisymmetric walls is presented in Fig.^. (a) Confinement into a film 
with antisymmetric boundary conditions enlarges the one phase region up to the prcwetting critical temperature. 
Since the wetting transition in binary polymer blends occurs far below the unmixing critical temperature in the bulk 
the effect is quite pronounced. The temperature region between the prewetting critical point and the triple point is 
about 11% of the wetting transition temperature. This value strongly depends on the details of the structure at the 
walls. The stronger the wetting transition the larger are the prewetting lines and the more extended is the region of 
the two miscibility gaps. The phase diagram of the bulk and a film with symmetric walls are displayed for comparison 
in FigJ^. The symmetric film has the same thickness as the antisymmetric film and the monomer-wall interactions at 
one wall are identical and attract the A-component. While the prewetting at the wall which prefers the A-componcnt 
leads to a two phase region in the antisymmetric case there is only a change in curvature of the binodal detectable in 
the symmetric case. 

Panel (b) of Fig.^j presents the phase diagram as a function of temperature and exchange chemical potential. In 
the antisymmetric case A/i cocx = up to the triple temperature. There, two coexistence lines emerge which are the 
thin film analogies of the prewetting lines at the two walls. Since the monomer-wall interactions are short ranged 
the prewetting line in the bulk and the coexistence curves in the film deviate from the bulk coexistence value linearly 
(up to logarithmic corrections). [^9| They end in two critical points. Though the system is strictly symmetric with 
respect to exchanging A x± B phase coexistence is not restricted to A/j, = 0, and the coexisting phases are not related 
by the symmetry of the Hamiltonian. The coexistence curve of the symmetric film is shown for comparison. The 
coexistence value of the chemical potential is shifted to values disfavoring the component attracted by both walls. 
There is a change in the temperature dependence of the coexistence curve close to the wetting transition temperature, 
but the coexistence curve stays far away from the prewetting line. If the two lines intersected there would also be 
a triple point in the symmetric case. [ pl|j5C| Since, the shift of the chemical potential A/i is roughly proportional to 
the inverse film thickness (Kelvin equation) we expect a triple point to occur only for much larger film thicknesses. 
This is in accord with self-consistent field calculations. [^l| The typical distance I between the interface and the wall 
at coexistence is of order D/2 in the antisymmetric case, while it is only of the order R g In D/R g in the symmetric 
case. Hence, smaller film thicknesses are sufficient to study the interaction between the interface and the wall, and 
antisymmetric boundary conditions are computationally more efficient to investigate the wetting behavior. 

V. SUMMARY AND DISCUSSION. 

We have studied the phase diagram of a symmetric polymer mixture in a thin film with antisymmetric boundary 
conditions via large scale Monte Carlo simulations. The walls interact with monomers via a short range potential; the 
one wall attracts the A component and repels the B component while the interaction at the opposite wall is exactly 
reversed. The salient features of the phase diagram and its dependence on the film thickness as obtained by our MC 
simulations are in accord with the results of mean field theory. |9[ |l7| , [l8|| Fluctuations, which are neglected in the mean 
field calculations, do not modify the qualitative phase behavior. However, they give rise to a rich crossover behavior 
between Ising critical behavior, tricritical behavior and their mean field counterparts. This has been elucidated by 
phenomenological considerations and is qualitatively consistent with our simulation results. 

Since the critical point of the thin binary polymer film occurs at much lower temperature than the unmixing 
transition in the bulk, "bulk-like" composition fluctuations are only of minor importance. The dominant fluctuations 
of the composition of the film arise from capillary waves at the interface between the j4-rich and i?-rich regions in 
the film. The interaction between the walls and the interface is rather small, because it is mediated via the distortion 
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of the interface profiles at the walls and the strength of the interaction decreases exponentially with the distance. 
Hence, the interface is only very weakly bound to the minimum of the effective interface potential. These large 
fluctuations give rise to rather pro nounced corrections to scaling in our systems of limited size. However, using the 
cumulant intersection method |35| and the matching of the order parameter distribution onto the predetermined 
universal scaling function, [p3 we give evidence for the 2D Ising universal character of the critical points. The same 
strategy has proven computationally very convenient to locate the tricritical point as a function of the film thickness. 
p6| This technique allows us to locate the critical points of the confined complex fluid mixture with an accuracy of a 
few percent. 

Interface fluctuation do not only impart 2D Ising critical behavior onto the critical points, but they are important 
in the whole temperature range. Monitoring the probability distribution of the laterally averaged interface position 
we extract the effective interface potential g(l). Its dependence on the lateral system size yields direct evidence for 
the renormalization of the interface potential by interface fluctuations. Interface fluctuations lead to a broadening of 
the minima in the interface potential, a shift of the minima towards the center of the film, and to a relative reduction 
of the free energy of the broader minimum. This leads to a systematic overestimation of the triple temperature by 
the mean field calculations. 

Moreover, our simulations indicate that packing effects in thin films result in corrections of the order 1/D to the 
density of the film or to the effective Flory-Huggins parameter. Such corrections are likely to mask completely 
the subtle thickness dependence of the triple temperature and the triple temperature predicted by the mean field 
calculations. For short range interactions between walls and monomers the predicted shifts decrease exponentially 
with the film thickness D. However, power-law dependencies are expected for the case of long range (i.e., van der 
Waals) interactions between walls and monomers. 

The gross features of the phase diagram as well as our simulation and analysis techniques are not restricted to binary 
polymer fluids but generally apply to binary liquid mixtures in confined geometries. Moreover, mean field calculations 
indicate that for small deviation from perfectly antisymmetric boundary conditions a qualitatively similar phase 
behavior emerges. The stronger the first order wetting transitions at the boundaries, the larger deviations from 
antisymmetry are permissible without alternating the topology of the phase diagram. Hence, a thin binary film on 
a substrate against air/vacuum, where the substrate energetically favors one component of the mixture while the 
other component has an affinity to the air surface, is an experimental realization of the boundary conditions discussed 
here. Our findings also imply that ultrathin enrichment layers at one surface are unstable in the temperature range 
T wc t < T < T c . Such effects have been observed experimentally |5l| in polymeric films, although for a liquid-vapor 
transition instead of a liquid-liquid demixing. However, recent experiments have observed the wetting transition in 
binary polymer blends, [j52T , j53| 
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TABLE I. Compilation of the boundaries of the different regimes in the vicinity of the tricritical point and the correlation 



lengths at the crossover. The latter quantity gives an estimate of the system size required to observe the crossover in the Monte 
Carlo simulations. 
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FIG. 1. (a) Illustration of the different regimes for a second order and tricritical transition: 2DTMF: mean 
field tricritical behavior, 2DMD: mean field critical behavior, 2DI: two dimensional Ising critical behavior, and 
2DT: two dimensional tricritical behavior. The inset shows the temperature dependence of the order parameter 
in for r = —0.4 as calculated within mean field theory (see Eq.(|)). For t c - 1 < 16r 2 /3 2DMF behavior is found, 
while 2DTMF behavior is observed at larger distances from the critical point, (b) Dependence of the critical 
temperature t c on the distance r from the tricritical point. The curves correspond to different values of AD as 
indicated in the key. Thick lines, which bracket the behavior, correspond to t c = 7r 2 /5 (valid for small r) and 
t c = 7r 2 /9 (valid in the limit AD — > oo). The inset presents the binodals at fixed strength b — 4.44 of the wetting 
transition of the individual surface and several values of AD as indicated in the key. For choice of parameter 
b/c — 4A4: > 3exp(— AD/2) (and, hence, r > 0) there are two critical points for all values of the film thickness. 

FIG. 2. Density of blocked lattice sites normalized by the bulk value as a function of the distance from the 
wall at e = 0.06 and e w = 0.16 for film thicknesses D = 24 and D = 48. Note the strong packing effects at the wall 
for z < 5. For these parameters an interface is stabilized at the center of the film. The position of the interface 
fluctuates in the interval R g ~7<z<D — R g (cf. Figj^) The inset presents the normalized density averaged 
over the layers 5-8. This region is marked by the arrow in the main panel. 

FIG. 3. (a) Probability distribution of the composition for system size D = 8 and L — 80. The inverse 
temperatures are indicated in the key. Histogram reweighting has been applied to extrapolate the data along the 
coexistence curve. The shape of the distribution function changes from single-peaked to bimodal, but there is 
no indication of a third peak at <j> = 1/2. (b) Cumulant ratio (m 2 ) / (\m\) 2 along the coexistence curve A/i = 
for film thickness D = 8 and various lateral extensions L as indicated in the key. In the finite size scaling limit 
L — > oo, t = (e — e c )/e c — > 0, Lt finite, the cumulant intersection should occur at the value (m 2 )/(|m|) 2 = 1.072, 
highlighted by the horizontal straight line. Our estimate of the critical temperature e c = 0.0520(5) is indicated 
by the double arrow. The inset shows the distribution function of the order parameter - scaled to unit norm and 
variance - at our estimate of the critical temperature and compares the MC results to the universal distribution 
of the 2D Ising universality class, (c) Same as (a) but for system size D = 12 and L = 48. Note that there is a 
broad range of e where the distribution has three peaks, unlike the Ising model. This indicates the vicinity of the 
tricritical point, (d) Cumulant ratio {m 2 } / (\m\) 2 for film thickness D = 12 and various lateral extensions L as 
indicated in the key. Our estimate of the critical temperature e c = 0.0589(10) is indicated by the double arrow. 
The inset shows the distribution function of the order parameter at our estimate of the critical temperature and 
compares the MC data to the universal distribution of the 2D Ising universality class. 

FIG. 4. (a) Probability distribution of the composition for various film thicknesses as indicated in the key. The 
lateral systems size is L = 96. We have adjusted the interaction strength e such that the central peak is a factor 
1.2 higher than the outer peaks. In accordance with convention we have scaled the distributions to unit norm and 
variance. Circles mark the universal distribution of 2D tricritical transition, (b) Temperature dependence of the 
cumulants for D = 14 and lateral system sizes as indicated in the key. The horizontal line marks the cumulant 
value of the universal tricritical distribution, (c) Probability distribution of the composition at e t ri = 0.06151(50) 
scaled to unit norm and variance. The universal 2D tricritical distribution (from Wilding and Nielaba p(|) is 
shown for comparison 

FIG. 5. (a) Temperature dependence of the cumulant (m 2 ) / {\rn\) 2 for D = 24 and various system sizes as 
indicated in the key. The arrow marks the critical temperature range e c = 0.061(1). (b) Probability distribution 
of the composition scaled to unit norm and variance at our estimate of the critical temperature e = 0.061. Thin 
lines denote the results of the Monte Carlo simulations. Histogram reweighting has been applied to extrapolate 
the data along the coexistence curve. Circles show the universal distribution of the 2D Ising universality class, 
(c) Same as (a) but for film thickness D = 48. The inverse critical temperature is e c = 0.0625(10). (d) Same as 
(b) but for film thickness D = 48 and e c = 0.0625. 

FIG. 6. (a) Dependence of the effective interface potential on the lateral system size L in a thin film of width 
D = 24 at e = 0.065 (b) same as (a) but for D = 48 and e — 0.069. The inset presents an enlarged view of 
the minimum close to the wall. The scale on the abscissa corresponds to the distance between the wall and the 
interface in units of R g . 
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FIG. 7. Free energy difference per unit area and ksT of the localized and delocalized state as a function of 
the lateral system size. The symbols represent the MC data, while the solid lines are calculated from the effective 
interface Hamiltonian. The temperature was chosen such that Ag — > for L — > oo. 

FIG. 8. Temperatures of the critical points and the triple point as a function of the film thickness D. Open 
symbols mark the results of the finite size scaling analysis. We have applied a correction factor (1 + 0.85/D) - 
to account for the film thickness dependence of the density at the center (filled symbols). Dashed lines are only 
guides to the eye. The arrow on the right hand side marks the value of the wetting transition temperature 
obtained independently via the Young equation, pi} 

FIG. 9. (a) Phase diagram of a binary polymer blend (N = 32). The upper curve shows the binodals in the 
infinite system; the middle one corresponds to a thin film of thickness D = 48 and symmetric boundary fields 
e w = 0.16, which both prefer species A. The lower curve corresponds to a thin film with antisymmetric surfaces. 
The arrow marks the location of the wetting transition. Full circles mark critical points; open circles/dashed 
line denote the triple point, (b) Coexistence curves in the (T, A/^)-plane. Circles mark critical points, and the 
diamond indicates the location of the wetting transition temperature. It is indistinguishable from the temperature 
of the triple point. 
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Fig.l : M.Miiller and K. Binder, Interface localisation — delocalisation in a symmetric ... 
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Fig. 2 : M.Miiller and K. Binder, Interface localisation — delocalisation in a symmetric ... 
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Fig. 4 : M.Miiller and K. Binder, Interface localisation — delocalisation in a symmetric ... 
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Fig. 7 : M.Miiller and K. Binder, Interface localisation — delocalisation in a symmetric ... 
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Fig. 8 : M.Miiller and K. Binder, Interface localisation — delocalisation in a symmetric . . . 
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Fig. 9 : M.Miiller and K. Binder, Interface localisation — delocalisation in a symmetric polymer 
blend . . . 
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